Imprints of Nuclear Symmetry Energy on Properties 
of Neutron Stars 



Bao-An Li^, Lie-Wen Chen^'^, Michael Gearheart^, Joshua Hooker^, 
Che Ming Ko'^, Plamen G. Krastev^'*, Wei-Kang Lin^, William G. 
Newton^, De-Hua Wen^'^, Chang Xu^'^ and Jun Xu'^ 

^Department of Physics and Astronomy, Texas A&M University-Commerce, Commerce, Texas 
75429-3011, USA 

■^Department of Physics, Shanghai Jiao Tong University, Shanghai 200240, China 
^Cyclotron Institute and Department of Physics and Astronomy, Texas A&M University, 
College Station, Texas 77843-3366, USA 

'^Department of Physics, San Diego State University, 5500 Campanile Drive, San Diego, CA 
92182-1233, U.S.A. 

^Department of Physics, South China University of Technology, Guangzhou 510641, China 
^Department of Physics, Nanjing University, Nanjing 210008, China 

E-mail: Bao-An_Li(9Tamu-Coimnerce . edu 

Abstract. Significant progress has been made in recent years in constraining the density 
dependence of nuclear symmetry energy using terrestrial nuclear laboratory data. Around 
and below the nuclear matter saturation density, the experimental constraints start to merge 
in a relatively narrow region. At supra-saturation densities, there are, however, still large 
uncertainties. After summarizing the latest experimental constraints on the density dependence 
of nuclear symmetry energy, we highlight a few recent studies examining imprints of nuclear 
symmetry energy on the binding energy, energy release during hadron-quark phase transitions 
as well as the w-mode frequency and damping time of gravitational wave emission of neutron 
stars. 



1. Introduction 

Properties of neutron-rich matter are at the heart of many fundamental questions in both nuclear 
physics and astrophysics [T]. They are currently being explored experimentally by using a 
wide variety of advanced new facilities, such as, Facility for Rare Isotope Beams (FRIB), x- 
ray satellites and gravitational wave detectors. Most critical to understanding experimental 
observations of various phenomena using these facilities is the Equation of State (EOS) of 
neutron-rich nucleonic matter. The latter can be written in terms of the energy per nucleon at 
density p as 

E{p, 5) = E{p, 6 = 0) + £;sym(p)<5' + 0{5^) (1) 

where 5 = {pn — Pp)/{pp + Pn) is the neutron-proton asymmetry and Esym{p) is the density- 
dependent nuclear symmetry energy. The latter is a vital ingredient in the theoretical description 
of neutron stars where the proton fraction is uniquely determined by the symmetry energy 
Esym{p) in E], and of the structure of neutron-rich nuclei and reactions involving them 
[H El El [71 [8]. It is currently the most uncertain part of the EOS of neutron-rich nuclear 



matter, especially at super-saturation densities. Very recently, attempts have been made to 
extract meaningful constraints on the neutron star EOS, facilitated by neutron star observational 
data with sufficient accuracy, and from a wide enough range of neutron star phenomena 
associated with different density regions within the star [9l \T0\ lllj . This signals a new era 
in which the quality and quantity of astronomical data make constraining the microphysical 
properties of neutron star matter from sub-nuclear to super-nuclear densities a realistic endeavor. 
At the same time, the nuclear matter EOS is being probed with unprecedented accuracy 
by a variety of terrestrial nuclear experiments, see, e.g., refs. [SJ [T21 [131 HH [151 [ISl [T7] . 
Combining the observational and experimental data, together with the astrophysical and 
nuclear modeling required to interpret them, has the potential to impose stringent constraints 
on the nuclear symmetry energy at both sub- and super-saturation densities, see, e.g., refs. 

2. The key physics behind nuclear symmetry energy 

Why is the nuclear symmetry energy so uncertain? To help understand this question, it is 
useful to first recall the relationship between the symmetry energy and the isospin-dependence 
of strong interaction at the mean-field level. According to the well-known Lane potential 
[23], the neutron/proton single-particle potential Un/p{p,k,6) can be well approximated by 
Un/pip,k,6) ^ Uo{p,k) ± Usym{p, k)5 where the Uo{p,k) and Usymip,k) are, respectively, the 
isoscalar and isovector (symmetry) potentials for nucleons with momentum k in nuclear matter 
of isospin asymmetry 5 at density p. It was shown first by Brueckner, Dabrowski and Haensel 
\25\ [26] using K-matrices within the Brueckner theory in the 60's-70's, and more recently by 
Xu et al. |27| using the Hugenholtz-Van Hove (HVH) theorem |28] that the nuclear symmetry 
energy can be explicitly expressed as 

Esym{p) = IfcF + ^Usym{P,kF), (2) 

where t{k) = h^k'^/2 m is the kinetic energy and kp = (37r^p/2)^/^ is the nucleon Fermi 
momentum in symmetric nuclear matter at density p. Moreover, it was shown recently using 
the HVH theorem that the slope of nuclear symmetry energy at an arbitrary density p can be 
written as [271 [2^] 

dE^p) 1 djt + Up) , i(p{t + Uo) 3 I, ,^9U^ 

L{P)=^P \p = g \kp-kF + -^ \kp-kF + l^U,yra{p,kF) + —Q^\kp-kF- 

(3) 

The Esym{p) and L{p) are correlated by the same underlying interaction. The above expressions 
for Esym{p) and L{p) in terms of the isoscalar and isovector single-particle potentials are 
particularly useful for extracting the symmetry energy and its density slope from terrestrial 
nuclear laboratory experiments. On one hand, in many reaction models, such as transport 
model simulations of nuclear reactions, the EOS only enters indirectly into the reaction dynamics 
and affects the final observables through the single-nucleon potential Un/p{p, k, 5). Except in 
situations where statistical equilibrium is established and thus many observables are directly 
related to the energy E{p, 6) after correcting for finite-size effects, what is being directly probed 
in nuclear reactions is the single-nucleon potential Un/p{p,k,6). On the other hand, in many 
structure models, such as shell models, the direct input is also the single nucleon potential. 
Thus, the isospin-momentum-density dependence of the single nucleon potential Un/p{p, k, 6) is 
the key in determining the EOS of neutron-rich nuclear matter. Comparing model calculations 
with experimental data allows the determination of the single nucleon potential which then 
determines directly the Esym{p) and L{p). While the density and momentum dependence of 
the isoscalar potential Uo{p,k) has been relatively well determined up to about 4 to 5 times 



the normal nuclear matter density po using nucleon global optical potentials from nucleon- 
nucleus scatterings as well as kaon production and nuclear collective flow in relativistic heavy- 
ion collisions [IS [13], our current knowledge about the symmetry potential Usym{p,k) is rather 
poor especially at high density and/or momenta [H [301 ED E21 [331 [31]- In fact, while some 
models predict decreasing symmetry potentials albeit at different rates, some others predict 
instead increasing ones with growing nucleon momentum. Thus, the high density/mometum 
behavior of the symmetry potential Usymip,k) is still very uncertain. Experimentally, there is 
some constraints on the symmetry potential only at normal density for low energy nucleons up 
to about 100 MeV obtained from nucleon-nucleus and (p, n) charge exchange reactions |27] . 

Of course, besides the uncertain density and momentum dependence of the symmetry 
potential, different approaches used in treating many-body problems in the various models also 
contribute to the divergence of the predicted symmetry energy especially at supra-saturation 
densities. Nonetheless, there are several key and commonly used physics ingredients that can 
affect the predicted Egymip) in all theories. In the interacting Fermi gas model, the isovector 
potential at kp can be written as [33 [36] 



in terms of the isosinglet (T=0) and isotriplet (T=l) nucleon-nucleon (NN) interactions Vroifij) 
and Vriirij), and the corresponding NN correlation functions f'^^{rij) and f'^^{rij), respectively. 
Needless to say, if there is no isospin dependence in both the NN interaction and correlation 
function, then the isovector potential Usymikp, p) vanishes. The -Esym(p) thus reflects the 
competition between the NN interaction strengths and correlation functions of the isosinglet and 
isotriplet channels. However, effects of several key ingredients, such as the tensor force, isospin- 
dependence of short-range nucleon-nucleon correlations and spin-isospin dependent three-body 
force, are not yet well understood, see, e.g., refs. (311 [33 [351 [32] . Consequently, the predicted 
symmetry energies diverge rather widely, especially at supra-saturation densities. 

3. Current status of experimental constraints on nuclear symmetry energy 

While significant progress has been made in constraining the £'sym(/o) as a result of the great 
efforts made by many people in both the nuclear physics and the astrophysics community, 
experimental constraints obtained so far are still not tight enough to limit stringently model 
predictions even around the saturation density pQ. Moreover, it has been broadly recognized 
that significant model dependence exists for most of the observables studied so far. For instance, 
near pQ the symmetry energy can be characterized by using the value of i?sym(po) and the slope 
parameter L{po). The £'sym(po) is known to be around 28 — 34 MeV mainly from analyzing 
nuclear masses within liquid-drop models |10l[lTl[32]. The exact value of -Bsym(po) extracted this 
way depends largely on the size and accuracy of nuclear mass data base as well as whether/what 
surface symmetry energy is used in the analysis. The latest study by Liu et al. [13] by 
analyzing the atomic masses indicates that the value of Esymipo) is between 29.4 and 32.8 
MeV while the value of L{pq) is between 52.9 and 78.7 MeV. The density slope L{pq) has 
also been extensively studied but remains much more uncertain. Shown in the left window of 
Fig. 1 are the L(po) values extracted from studying several phenomena/observables in nuclear 
structures and reactions [271[Mlia[l6l[lZllMll9l[Mll5ll[Ml[531[Ml[55]- The values of 
L{pq) from these studies scatter between about 20 to 115 MeV while each individual study 
may have given a smaller uncertain range. Although it is very encouraging to note that the 
constraints seem to start converging around L(po) ~ 60 MeV and -EsymCpo) ~ 31 MeV, it is 
obviously important to further narrow them down with new measurements and model analyses. 
It is thus encouraging to see the most recent determination of Egymipo) = 30.5 ± 3 MeV and 
L{po) = 52.5 lb 20 MeV [56] from combining the constraints on L{po) and Esym{po) based on 





Figure 1. Left Window: Density slope versus tlie magnitude of the symmetry energy at normal 
density extracted from isospin diffusion and neutron/proton ratio of pre-equilibrium nucleon emissions 
within the Improved Molecular Dynamics (ImQMD-2009) model [3111131135], isospin diffusion within the 
Isospin-Dependent Boltzmann-Uehling-Ulenbeck (IBUU04- 2005) [3Z1I3H]; isoscaling [32], energy shift of 
isobaric analogue states within liquid drop model (IAS+LDM-2009) [SD], neutron-skins of several heavy 
nuclei using the droplet model (DM) or the Skyrme-Hartree-Fock (SHF) approach [511 [52l [53] , pygmy 
dipole resonances (PDR) in ^"^Pb, ^®Ni and ^'^^Sn [M] [55], and the nucleon global optical potentials 
(GOP) [27]. Taken from ref. [27]. Right Window: Comparisons of the scaled symmetry energy as 
a function of the scaled total density n/riQ for different approaches and the experiment. Left panel: 
Symmetry energies from the MDI parametrization of Chen et al. [471 165] for T = and different 
L{pq), controlled by the parameter x (dotted, dot-dashed and dashed (black) lines); from the quantum- 
statistical (QS) approach including light clusters for temperature T = 1 MeV (solid (green) line), and 
from the relativistic mean-field (RMF) model at T = including heavy clusters (long-dashed (orange) 
line). Right panel: The internal scaled symmetry energy in an expanded low-density region. Shown are 
again the MDI curves and the QS results for T = 1,4, and 8 MeV compared to the experimental data 
using the entropy from the nuclear statistical equilibrium (NSE) model (solid circles) and the results of 
the self-consistent calculation (open circles). Taken from ref. [55] . 



recently extracted EsymiPA = 0.1 fm~^) [IS] [551 ES] with those from recent analyses [53] of 
neutron skin thickness of Sn isotopes in the same SHF approach. At very low densities, the 
formation of clusters is energetically favored. The symmetry energy of the clustered matter is 
expected to be different from that of uniform matter, see, e.g., refs. [601 [611 [62] . Indeed, it was 
recently found in experiments by Natowitz et al. that the symmetry energy of clustered matter 
remains finite as the average density approaches zero [63]. As shown in the right window of 
Fig. 1, conventional theoretical calculations of the symmetry energy of uniform matter based 
on mean-field approaches (e.g., the MDI interaction [64[ Wl\ I65j ) fail to give the correct low- 
temperature, low-density limit that is governed by correlations, in particular by the appearance 
of bound states. A recently developed quantum statistical (QS) approach [62] that takes the 
formation of clusters into account predicts symmetry energies that are in very good agreement 
with the experimental data. At supra-saturation densities, the situation is much worse. In 
fact, even the trend of the symmetry energy, namely, whether it increases or decreases with 
increasing density, is still controversial partially because of the very limited data available and 
the less known clean probes compared to the situation at sub-saturation densities. Among all 
observables studied so far, the neutron-proton differential flow [66) . 7r~ /tt~^ ratio [67| and the 
squeeze-out of the neutron/proton ratio ^8] especially at high transverse momenta in heavy- 
ion reactions are among the most promising probes of the high density behavior of the nuclear 
symmetry energy. Unfortunately, the conclusions are still model dependent [69| \70 \ \T\ ] [72] [73]. 



4. Impacts of nuclear symmetry energy on properties of neutron stars 

It is well known that many critical issues in astrophysics, such as the composition and the 
cooling rate of proto-neutron stars as well as the core-crust transition density, properties of the 
pasta phase, the mass-radius correlation and the moment of inertia of neutron stars, all depend 
sensitively on the Ssym(p)- Many research papers and review articles have been written on these 
topics, see, e.g., refs. [2l[3l[71]. Most of these studies have concentrated on astrophysical effects 
of the symmetry energy at sub-saturation densities. However, the nuclear symmetry energy 
at supra-saturation densities has long been considered by some astrophysicists as the most 
uncertain one among all properties of dense neutron-rich nucleonic matter |75 [ I76 [ [77]. More 
efforts are thus need to explore many possible astrophysical effects of the symmetry energy at 
supra-saturation densities [75]. In this section, as examples we highlight a few recent studies 
on effects of the symmetry energy, especially its high-density behavior, on the binding energy 
of neutron stars, total energy release during mini-collapse triggered by hadron-quark phase 
transitions inside neutron stars, and the w-mode of gravitational waves from neutron stars. 

4-1- Effects of symmetry energy on the binding energy of neutron stars 




Figure 2. Left Window: Correlation between the slope parameter L{pq) of the symmetry energy and 
the compactness parameter (3 — GMq/Rc^ for the sample of EoSs. Right Window; The fractional 
gravitational binding energy of the neutron star BE/Mq versus the compactness parameter /3 for the 
sample of EoSs. The dashed line indicates the analytic Newtonian result. Taken from ref. [23]. 

The binding energy of a neutron star is given by [79] BE = Mq — Mbarj where Mq is the 
gravitational mass of the neutron star measured from infinity [80], i.e., 

Mg = I 4Trr'^e{r)dr (5) 
Jo 

with e(r) being the energy density at radius r. The Mbar is the baryonic mass of the neutron 
star, namely, the mass when all the matter in the neutron star is dispersed to infinity |79j . It 
can be calculated from Mbar = NMb, where Mb = 1.66 x 10~^^g is the nucleon mass and is 
the total number of baryons [79], i.e., 

-1/2 

PB{r)dr (6) 



N 



A-nr'^ 



2m[r) 



with pb{i~) being the baryon density profile of the neutron star. For a given EOS, the binding 
energy can be readily obtained by first solving the Tolman-Opperheimer-Volkov (TOV) equation. 
Newton et al. [23] examined the correlation among the symmetry energy slope L{po), pressure 
at normal density P{nQ), compactness of neutron stars /? = GMq/R(? and the binding energy. 
As shown in Figs. El for a wide range of EoSs, the binding energy correlates roughly linearly with 
the compactness of the neutron star /3, which in turn correlates strongly with the slope of the 
nuclear symmetry energy L{pq) predicted by the EoSs. Therefore one expects to see a correlation 
between L{po) and the binding energy of the neutron star, and consequently the baryon mass, 
given that the gravitational mass is fixed. Thus, laboratory constraints on L can be used to 
extract the constraint on the baryon mass of neutron stars. This approach was used to constrain 
the lower mass member of the double pulsar binary system, PSR J0737-3039B. Comparing with 
independent constraints derived from modeling the progenitor star of J0737-3039B up to and 
through its collapse under the assumption that it was formed in an electron capture supernova, 
it was found that the two sets of constraints are consistent only if L(po) < 70 MeV [23]. 




Figure 3. Left Window: Density dependence of nuclear symmetry energy from the MDI 
interaction with parameter x = and —1. Right Windows: EOSs of pure npefi matter 
(nucleonic), hyperonic matter (MDI Hyp-R interaction) and hybrid stars for the bag constant 
5^/4 = 180 MeV and 170 MeV with x = (a) and x = -1 (b). 



The total energy release Eg during the hadron-quark phase transition is the difference 
of binding energy before and after the micro-collapse. It reduces to the difference in the 
gravitational masses for the hadronic (h) and hybrid (q) configurations of the neutron star 
as a result of baryon number conservation, namely. Eg = Mq ^ — M^q- Thus, the Eg should 
depend on the EOSs of dense matter either in the hadronic or quark phase. For the neutron-rich 
hadronic matter, the most uncertain part of its EOS is the density dependence of the nuclear 
symmetry energy, while for the quark matter within the MIT bag model, the most uncertain 
part of its EOS is the bag constant. The relative effects of the density dependence of the nuclear 
symmetry energy and the bag constant on the total energy release due to the hadron-quark phase 
transition in NSs were recently examined by Lin et al. [81]. Shown in Fig. [3] are several model 
EOSs for hybrid stars obtained from an isospin- and momentum-dependent effective interaction 
(MDI) [Ml [82] for the baryon octet, the MIT bag model for the quark matter [831 Hi] and the 
Gibbs construction for the hadron-quark phase transition [851 186j. For comparisons, the pure 
npep (labeled as nucleonic) and hyperonic (labeled as MDI Hyp-R interaction) EOSs are also 
included. As it is well known, the appearance of hyperons and the hadron-quark phase transition 
softens the EOS of neutron star matter. The EOSs with x = and x = —1 are about the same 



below and around the saturation density. However, it is interesting to see that the starting point 
and the degree of softening due to the appearance of hyperons are sensitive to the Egymip) at 
high densities. Moreover, the i?sym(/') also affects appreciably the starting point of the hadron- 
quark mixed phase, especially with the larger bag constant. Nevertheless, it is obvious that 
the starting point is much more sensitive to the bag constant B for a given symmetry energy 
parameter x. These features are consistent with those first noticed by Kutschera et al. [87j . 
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Figure 4. Left Window: Total energy release due to the hadron-quark phase transition as a function 
of the baryonic mass of a neutron star. Right Window: Damping time versus frequency of the first 
w-mode for hybrid stars. 



Shown in the left window of Fig. U] is the energy release as a function of the baryonic mass 
-^bar/-^© of a neutron star. It is seen that the energy release increases with Mbar/-^© and 
is higher with a smaller (B^^'^ = 170 MeV) bag constant B but a stiffer {x = —1) symmetry 
energy. Effects of varying the bag constant B are obviously more significant than varying the 
symmetry energy parameter x, especially on the core mass and thus the energy release. The 
variation of the bag constant B also affects appreciably the radii of hybrid stars. On the other 
hand, the variation of the symmetry energy parameter x only has appreciable effects on the radii 
of both hyperonic and hybrid stars. It's effects on the energy release is much smaller than the 
bag constant B. 

4-3. Imprints of symmetry energy on the w-mode of gravitation waves 

The search for gravitational wave signals is another forefront in astrophysics and cosmology. 
Gravitational waves are tiny disturbances in space-time and are a fundamental, although not yet 
directly confirmed, prediction of General Relativity. They would open up an entirely new non- 
electromagnetic window for probing the physics that is hidden or dark to current electromagnetic 
observations. There are many possible sources of gravitational waves. For instance, the hadron- 
quark phase transition can trigger various oscillation modes of neutron stars. When coupled 
with rotation, these oscillations generate gravitational waves. The energy release discussed 
in the previous subsection determines the maximum strain amplitude of emitted gravitational 
waves. Moreover, the strain amplitude of gravitational waves from pulsars depends on their 
quadrupole deformation which, in turn, is determined by the EOS of neutron-rich nucleonic 
matter. Furthermore, according to Chandrasekhar & Ferrari [88], the frequency and damping 
time of the u;-mode of gravitational waves from oscillating neutron stars are described by a 
unique second-order differential equation with the EOS as the only input. Indeed, it was recently 
shown that the strain amplitude from elliptically deformed pulsars as well as the frequency and 
damping time of the tt;-mode of oscillating neutron stars are all strongly influenced by the density 



dependence of the nuclear symmetry energy [891 [901 [9T1 [921 [93] . Thus, an accurate determination 
of the nuclear symmetry energy is also important for detecting gravitational waves using various 
ground-based and satellite-based detectors, such as, LIGO [93] and LISA [95]. As a quantitative 
example, shown in the right window of Fig. [J] is the frequency versus damping time of the first 
axial w-mode. For the softer EOS with S^/^ = 170 MeV, the effect of the symmetry energy is 
small. The density dependence of the symmetry energy can, however, significantly affect the 
frequency of the ■w-mode for the stiffer EOS with B^^'^ = 180 MeV as a result of the higher 
hadron-quark transition density. In this case, the symmetry energy in the hadronic phase then 
also plays an important role in determining the structure of NSs. Nevertheless, comparing the 
frequencies of the w-mode for the same x parameter but different values of B, it is clear that 
the bag constant has a much stronger effect as it changes the underlying EOS of dense matter 
more significantly. 

5. Summary 

In summary, important progress has been made in recent years in constraining the symmetry 
energy, especially around and below the saturation density, with terrestrial nuclear laboratory 
data. Nevertheless, the field is still at its beginning. While a number of potentially useful probes 
of the symmetry energy, especially at supra-saturation densities, have been proposed, available 
experimental data are mostly for reactions with stable beams. Coming experiments with more 
neutron-rich nuclei at several advanced high energy radioactive beam facilities are expected to 
improve the situation dramatically. Thus, more exciting times are yet to come. Implications 
of the partially constrained symmetry energy on some astrophysical phenomena have also been 
explored. On the other hand, impressive progress has also been made in extracting information 
about the EOS from astrophysical observations. Hopefully, working together and combining 
results from both fields will allow us to finally pin down the EOS of neutron-rich nuclear matter 
at both sub- and super-saturation densities in the near future. 
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